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Abstract 

The mirror system of the HERA-B RICH consists of two spherical and two planar 
mirrors, composed of altogether 116 mirror segments. Analysis of displacements of 
the Cerenkov ring center relative to the charged particle track, for given spherical- 
planar segment pairs, leads to accm'ate information regarding the orientation of 
individual mirror segments. The method is described and the effect of applying 
the required corrections on the Cerenkov angle resolution of the HERA-B RICH is 
discussed. 
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1 Introduction 



HERA-B (l) was a fixed target experiment (Fig. 1) at the HERA storage ring 
at DESY in Hamburg. The experiment used 920 GeV protons from the beam 
halo and a set of eight thin ribbons, of different materials, as targets. The in- 
teraction rate was adjusted by moving the targets in or out of the beam halo 
(12). The experiment utilized a forward spectrometer capable of measuring in- 
teraction rates up to 40 MHz. The spectrometer consisted of a dipole magnet, 
a vertex detector (jSl) upstream and a main tracking system downstream of 
the magnet (|J; |5|). Particle identification was performed by a Ring Imaging 
Cerenkov (RICH) detector (6), an electromagnetic calorimeter (j) and a muon 
detector system (jsl). In addition, the experiment included a sophisticated hard- 
ware trigger for lepton track pairs to record leptonic decays of J/ip particles. 
The large acceptance of the spectrometer coupled with high-granularity par- 
ticle identification devices and a precision vertex detector allowed for detailed 



studies of multi-particle final states (l9l)-(ll6|). By using targets of different ma- 



terials, HERA-B was also able to study the dependence of various properties 
of proton-nucleus interactions as a function of atomic number. 



The identification of pions, kaons and protons was performed by the RICH de- 
tector (p). The HERA-B RICH used atmospheric pressure C4F10 as Cerenkov 
radiator (n=l. 00137). The focusing of Cerenkov light was achieved with two 
spherical mirrors, tilted by 9° in opposite directions (Fig. 1). Two planar mir- 
rors then reflected the light to photon detectors at the top and bottom of the 
vessel containing the radiator gas. For the detection of Cerenkov photons multi 
anode PMTs (Hamamatsu R5900) were used. The inner part of the photon 
detector surface was equipped with 16-channel PMTs (16 x 4 mm x 4 mm) 
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and the outer region with coarser granularity had 4-channel PMTs (4x8 mm 
X 8 mm). To overcome the loss of photons due to inactive space between PMT 



photocathodes, a demagnifying lens system (1171 ) was placed in front of each 
PMT. 
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Fig. 1. A side view of the HERA-B detector. Photon paths in the RICH counter 
are indicated. 



The particle identification capabilities of a RICH counter are determined by 
the resolution of the measured Cerenkov angle, which is given by the two main 
parameters of a RICH counter, the Cerenkov angle resolution due to a single 
photon and the number of detected photons per Cerenkov ring. The measured 
average number of detected photons for particles approaching the speed of 
light amounts to 33. It is in good agreement with the value expected from the 
data available on the quantum efficiency, mirror reflectivity, and transmissions 
of the vessel window and of the optical system (6). 

The main contributions to the single photon resolution (r.m.s.) come from the 
photon detector granularity (0.50 mrad and 0.93 mrad for the regions covered 
by finer and coarser granularity PMTs respectively) and the dispersion in the 
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radiator medium (0.33 mrad). The optical error (0.25 mrad) includes contri- 
butions from spherical aberration, mirror quality, and mirror alignment. The 
contribution of multiple scattering in the RICH counter ( p^(G^/c) ) becomes 
important at low momenta. The resulting expected single photon resolution, 
0.65 mrad © and 1.02 mrad © ^^nT/'^^ for the regions covered by the 

p (GeV/c) p (GeV/c) 

two types of PMTs, does not include the contribution from the uncertainty 
in the track direction, which is given by other components of the HERA-B 
detector. 

In order to reach and maintain the optimal performance of the RICH counter, 
elaborate alignment and calibration methods have to be used. For the optical 
system of the HERA-B RICH, a calibration method was used which is based 



on a procedure originally developed on simulated data (llSl ). 



2 Calibration of the optical system 

The main imaging device of the HERA-B RICH is a spherical mirror placed 
inside the radiator vessel with the center of the sphere near the target and a 
radius of curvature of 11.4 m. The mirror, a 6 m by 4 m rectangular cutout 
from the sphere, consists of 80 full or partial hexagons (see Fig. 2). To achieve 
a focal surface outside of the main particle flux (±160 mrad vertically), the 
mirror is split horizontally, and both halves are tilted by 9° away from the 
beam-line. A set of two planar mirrors, composed of 18 rectangular elements 
each, translates the focal surface to the photon detector area above and below 
the radiator vessel (see Fig. 1). The 116 mirror segments are mounted on rigid, 
low mass support structures inside the radiator volume and can be individually 
adjusted by stepper motors from the outside. 
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6 m ► 

Fig. 2. Distribution of spherical mirror polygons. The holes in the array are for the 
proton and electron beam pipes. 



All mirrors were first aligned after installation by surveying them inside the 
vessel. During the data taking periods, the mirror system was calibrated by 
making use of recorded events. By comparing the charged particle track di- 
rection, obtained from the Cerenkov rings due to a particular spherical-planar 
mirror pair, to the track direction obtained from other detectors of HERA-B, 
the calibration parameters of individual mirror segments, as well as those of 
the entire RICH counter, could be extracted. For this purpose, various data 
sets have been used. With the magnetic field turned off, the direction of the 
straight tracks was accurately given by the target wire position and the cen- 
troid of the cluster in the electromagnetic calorimeter. With the magnetic field 
turned on, the tracks were determined by the tracking system. In order to re- 
duce the uncertainty in track direction due to multiple Coulomb scattering, 
only those tracks belonging to particles with energy above 5 GeV were used. 



5 



2.1 The calibration method 



Assume that one or both mirrors in a particular spherical-planar mirror pair 
are not well aligned. In such the measured Cerenkov rine;, due to 

photons reflected on that pair, will be displaced relative to the direction of 
the charged particle, which is taken to be reflected on ideal mirror positions. 
For small displacements a, the azimuthal dependence of Cerenkov angle for 
photons on a given ring is parameterized as (Fig. 3): 

dc = + a cos(0c — 0o) = do + A$ cos 0c + sin (pc, (1) 



where is the nominal value of the Cerenkov angle. The parameters A$ = 
a cos 00 and AA = asin0o roughly correspond to rotations of the mirrors 
around vertical and horizontal axes. A rotation of a spherical mirror segment 
by 6 around the vertical and horizontal axes results in A$ = 6 and AA = 6, 
respectively, while the same rotation of the planar mirror results in A$ ^6/2 
and AX ^ 6/2 (see Fig. 4). 









\ J' ^ 




Fig. 3. For misaligned mirrors, the center of the measured ring C is displaced 
relative to the extrapolated charged particle track direction C (left). In such a case, 
the measured Cerenkov angle 9c depends on the azimuth 0c of the photon hit (right). 
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The method obviously rehes on accumulating a sufficient number of photons, 
which have been reflected on a particular spherical-planar mirror pair. For 
each track-photon pair, the photon has been traced from two points on the 
charged particle track to the photon hit position (Fig. 4). The flrst point 
is the particle entry into the radiator, the second is directly in front of the 
spherical mirror. If both rays at a given azimuthal angle are reflected from the 
same spherical-planar mirror pair, such a photon hit is taken as a valid data 
point on the Cerenkov ring, relevant for that mirror pair. In other words, a 
photon hit is valid for our analysis if the corresponding photon would have 
been reflected from the same spherical-planar mirror pair, regardless of the 
point on the charged particle trajectory from which it might have originated. 
In the following we refer to such hits as calibration hits. 





Planar mirror 

^ ' Spherical mirror 




E track Ja ^ 
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Fig. 4. For calibration of mirror orientations, only those photon hits have been 
used, for which the corresponding photons could have been reflected only from a 
given spherical-planar mirror pair, regardless of the point on the particle track, from 
which the photon was emitted. 
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Figure 5 shows the distribution of such hits for two particular combinations 
of a spherical and two planar mirrors. The two-dimensional histograms on 
the left hand side of the figure represent the number of calibration hits as a 
function of Cerenkov angle difference A^c = dc — dc{p) ^'^d azimuthal angle 0c- 
The difference to the nominal Cerenkov angle {p) of a pion at the measured 
momentum p is chosen in order to suppress the momentum dependence; the 
pion hypothesis is chosen since the majority of tracks correspond to pions. An 
accumulation of hits at about A6c = can be seen, with some modulation 
dependent on azimuthal angle 0c- For each of 50 slices in (pc, the distribution 
was fitted with a Gaussian for the peak and a polynomial background (Fig. 6). 
The Cerenkov peak position as a function of azimuthal angle (right hand side 
of Fig. 5) is then fitted with the function A$ cos 0c + AA sin 0c + C, and the 
rotation angles A$ and AA are obtained for the given mirror paij^. 

Close to mirror boundaries, parts of a ring could be shared by different com- 
binations of spherical and planar mirrors, which permits relative calibration 
of adjacent mirrors. 

2.2 Extraction of calibration parameters of individual mirror segments 

The measured displacement (A$, AA) for a pair of mirrors is equal to the sum 
of unknown contributions from the corresponding spherical and planar mirror 
segments k and j: 

A^kj = A^f^ + A$f , (2) 

^ The parameter C is a constant to account for a possible bias in the Cerenkov 
angle measurement. 
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Fig. 5. Distribution of hits in the A^c) <Pc plane for two combinations of a spherical 
and two planar mirror segments. Raw data are shown on the left side, and the ring 
peak position from the fit in each (f)c slice is shown on the right. 

AAfej = AAf'^ + AA^'. (3) 



This represents a system of 2m linear equations for the 2n unknown parame- 
ters, where m is the number of all measured segment combinations {k,j), and 
n is the number of mirror segments. As any given mirror segment contributes 
to more than one combination, there are more equations than there are un- 
knowns. Arranging the measured A$fcj and AA^j into vectors bf and bf, with 
dimensions equal to m, the system of equations may be written for each of 
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Fig. 6. Distribution of hits in A^c for one of the shces. The result of the fit is 

superimposed. 

the detector halves as 

n n 

bf = J2Aauf, bf^J2^ii^t, i^l,...,m, (4) 

1=1 1=1 

where uf and are the unknown contributions to the rotation (A$, AA) of 
a particular mirror segment. Note that the matrix An is of a particular simple 
form; in a given row it only has 1 at two places (given by the indices of the 
spherical and planar segments, Eqs. 2 and 3), and is otherwise equal to zero. 

Note that we have to add an additional condition to fully determine the sys- 
tem: if all spherical mirrors are turned by S, and all planar mirrors by ^ —25, 
one arrives at the same ring displacements (vectors bf and bf). Rather than 
fixing the displacement of one mirror in the lower and one in the upper half 

to zero (one of the possibilities), we require of the solution to minimize the 
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necessary mirror readjustments, 

E<-E< = o, E<-E< = o. (5) 

The first sum in each of the equations runs over all spherical mirror segments 
and the second over all planar ones in a given detector half. We include these 
two equations in the system of equations 4, and get two systems with m + 1 
equations, 

n n 

bf = Y.A',,uf, bf = Y.A',,ut, 2 = l,...,m + l, (6) 

1=1 1=1 

where A' replaces A to account for the additional equation. 

The systems of equations, Eq. 6, are solved by requiring that the properly 
weighted sum of squares of deviations of the left hand side from the right 
hand side for each of the two systems is minimal, 

2 = i^) 

1=1 

where the terms are weighted by the inverse square of the error cxj in the 
measuremenjf] of 6j. Here ui and bi denote either uf and bf or and b^. The 
resulting linear system, 

n m+1 A> Al m+1 A' 

E E = E -fb. (8) 

1=1 i=l "i i=l "i 



is readily solved, 

n m+1 Al 

Uk = E -4b^■ (9) 

j=l i=l 



^ For the two additional equations 5, a value of ui = 0.01 mrad was assumed; no 
influence was found when this value was varied. 
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Here we have defined a new matrix 

m+l Af Al 

Bij = E (10) 

i=l '^i 



which is symmetric, so that its inverse B -"^ is easy to calculate. The resulting 
errors on Uk are given by 



m+l n A I 

t — T ,t — T , _ 1 , 



< = E(E(5-%^f- (11) 



1=1 j=i 



2.3 Calibration results 



Having solved the system of equations for angular displacements of individual 
mirror segments, the quality of the new alignment is checked on the data, by 
examining the Cerenkov angle resolution before and after applying the new 
cahbration parameters, A$^^'^, AA^^'^, A$^\ and AA^'. The result is shown in 
Fig. 7; the improvement in resolution is clearly visible. 

To check the resulting alignment of the optical system, we have reanalyzed 
the data after having applied the calculated corrections. As expected, the new 
corrections were consistent with zero. 

To study possible systematic effects, we have investigated different data sets, 
recorded under different conditions. We found good agreement of the results 
with two different gas radiators, when freon (6*0=52 mrad) was replaced with 
nitrogen {9c— 25 mrad). A similar analysis was also performed on the data 
recorded without magnetic field. Again the analysis yielded results that are 
consistent with the values deduced from the data with magnetic field. We also 
found that the alignment parameters did not change over extended periods of 
time, which testifies to stable mirror positions. 
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Fig. 7. The distribution of photon hits with respect to the corresponding Cerenkov 
angle difference A^c before (a) and after (b) applying the correction obtained from 
the calibration of mirrors. 

Possible systematic effects were also checked by using Monte Carlo generated 
events, where all mirror segments were assumed to be perfectly aligned. Al- 
though the resulting parameters are all consistent with zero, small systematic 
effects at the level of 0.1 mrad could not be excluded due to hmited statistics. 



3 Conclusions 

A method was developed for determination of the ahgnment of the RICH 
counter relative to other parts of the spectrometer, either the system of track- 
ing chambers or the electromagnetic calorimeter. The method was tested on 
various sets of real data, recorded with and without magnetic field, as well as 
on simulated data. By applying alignment corrections for each mirror segment 
as derived by this method, a significant improvement in the resolution of the 
Cerenkov angle measurement could be obtained. 

The HERA-B experiment finished data taking in spring 2003. During its five 
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years of operation, the HERA-B RICH has proved rehabihty and stabihty of 
all its components and especially of the multi anode PMTs. The low noise, high 
rate capability and excellent long term stability of these devices enabled excel- 
lent operation in the hostile environment of a hadron machine. The HERA-B 
RICH identifies pions, kaons and protons essentially in the entire kinematic 
range of the HERA-B experiment with the identification efficiencies as large as 
90% and mis-identification probability at the 1% level (jsjIlO). With kaon and 
proton identification the combinatorial background is in some cases reduced 
by more than 3 orders of magnitude. Several physics analyses would not have 
been possible without the excellent performance of the RICH particle identi- 



fication system (|19 
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